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We systematically investigate the ground-state and the spectral properties of antiferromagnets 
on a kagome lattice with several common types of the planar anisotropy: XXZ, single-ion, and 
out-of-plane Dzyaloshinskii-Moriya. Our main focus is on the role of nonlinear, anharmonic terms, 
which are responsible for the quantum order-by-disorder effect and for the corresponding selection 
of the ground-state spin structure in many of these models. The XXZ and the single-ion anisotropy 
models exhibit a quantum phase transition between the q = 0 and the a/3 Xa/3 states as a function 
of the anisotropy parameter, offering a rare example of the quantum order-by-disorder fluctuations 
favoring a ground state which is different from the one selected by thermal fluctuations. The 
nonlinear terms are also shown to be crucial for a very strong near-resonant decay phenomenon 
leading to the quasiparticle breakdown in the kagome-lattice antiferromagnets whose spectra are 
featuring flat or weakly dispersive modes. The effect is shown to persist even in the limit of large 
spin values and should be common to other frustrated magnets with flat branches of excitations. 
Model calculations of the spectrum of the S' = 5/2 Fe-jarosite with Dzyaloshinskii-Moriya anisotropy 
provide a convincing and detailed characterization of the proposed scenario. 

PACS numbers: 75.10.Jm, 75.30.Ds, 75.50.Ee, 78.70.Nx 


I. INTRODUCTION 

Kagome-lattice antifferomagnets are iconic in the field 
of frustrated magnets, comprising a number of model 
systems whose classical ground states are massively de¬ 
generate, giving rise to an extreme sensitivity to sub¬ 
tle symmetry breaking effects,^!^^ to fractional magneti¬ 
zation plateaus ,1^ to a strongly amplified role of sec¬ 
ondary inter act ions ,1^ and to an emergent hierarchy of 
energy scales in their dynamics.A crucial role of the 
non-linear, anharmonic terms in the so-called order-by¬ 
disorder ground-state selection by thermaP or quantum 
fluctuation^^ in the kagome-lattice antifferomagnets has 
been recognized for some time. Recently, an accurate, 
systematic treatment of the quantum order-by-disorder 
effect in the anisotropic versio ns of the model has re¬ 
ceived a significant developmentHowever, much less 
is known about the role of such terms in the excitation 
spectra of frustrated magnets and only recently a rather 
dramatic picture has begun to emerge.tm 

Usually, the nonlinear terms in antiferromagnets are 
necessary to describe interactions of magnons and, while 
their role is, ge nerally, more significant when frustration 
is presentthey still lead to effects that are relatively 
minor in the large-S' limit, i.e., constitute a 1/S con¬ 
tribution compared with the classical energy scale JS. 
However, in a wide class of highly-frustrated systems, 
including the kagome-lattice antiferromagnets, the non¬ 
linear anharmonic terms are responsible for the phenom¬ 
ena that are much more dramatic. 

First effect concerns systems in which contenders for 
the ground state form a highly degenerate, extensive 
manifold of states, and neither the classical energetics 
nor the harmonic fluctuations are able to select a unique 


ground state from that manifold. h l ^ l ^ l In that case, the 
selection role is passed onto the nonlinear terms pro¬ 
viding a variant of the quantum order by disorder ef¬ 
fect. Such is the case of the three models considered in 
this work, Heisenberg, XXZ^ said single-ion anisotropy 
models, while in the case of the Dzyaloshinskii-Moriya 
(DM) anisotropy, a unique state is selected already on 
the classical leveL^iS 

The second effect is less-studied, but is equally strik¬ 
ing. We demonstrate, that the nonlinear terms lead to 
spectacularly strong quantum effects in the dynamical 
response of the flat-band frustrated magnets, even in 
the ones that are assumed nearly classical.!^ The resul¬ 
tant spectral features invoke parallels with the quasipar¬ 
ticle b reakd own signatures in quantum spin- and Bose- 
liquids,!!^^ which exhibit termination points and broad 
continua where single-particle excitations are no longer 
well-defined. In the present case, the origin of such fea¬ 
tures is in the near-resonant decay into pairs of the flat 
modes, facilitated by the nonlinear couplings. The effect 
is strongly amplified by the density of states of the flat 
modes and can be shown to persist even in the large-5' 
limit, defying the usual 1/5 suppression trend and chal¬ 
lenging a conventional wisdom that such drastic phenom¬ 
ena can only occur in an inherently quantum system. 

In the pr esent study, we expand the analysis of our pre¬ 
vious work^^^ and offer an expose of the 1/5 formalism 
for several common anisotropic extensions of the nearest- 
neighbor Heisenberg model on the kagome lattice that are 
also relevant to real materials. The presented approach 
to the nonlinear spin-wave theory can be applicable to 
the other, more complicated forms of the kagome-lattice 
Hamiltonians as well as to a broader class of frustrated 
spin systems on the non-Bravais lattices. We also provide 
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a useful extension of our approach to the perturbative 
treatment of small perturbations to the main Hamilto¬ 
nian, such as the next-nearest superexchanges J 2 . 

For the ground-state consideration, we demonstrate a 
quantum phase transition between the q = 0 and the 
\/3 X \/3 states as a function of anisotropy parameter in 
two models: XXZ model, studied previously,!^ and the 
single-ion anisotropy model. While the effect is similar in 
both models, the energy scale associated with it is shown 
to be different, in agreement with the understanding 
that the degeneracy-lifting interaction is associated with 
a high-order topologically-nontrivial spin-flip processes, 
which are different in the two models. Nevertheless, both 
model cases present a rare example of the ground-state 
selection that is different from the choice of the thermal 
fluctuations, which favor the \/3 x \/3 struct ure for any 
value of the anisotropy parameter 

For the spectral properties of the kagome-lattice an- 
tiferromagnets, we offer a detailed consideration of the 
decay-induced effects in the DM anisotropy model with 
J 2 , the mod el that closely describes the S = 5/2 Fe- 
jarosite.l^^^ We also present a general analysis of the 
near-resonant decays in the flat-band frustrated antifer- 
romagnets, which suggests that the dramatic modifica¬ 
tions in the spectrum due to this phenomenon must per¬ 
sist in the large-S' limit. While the core of our presen¬ 
tation is aimed at a common and realistic extension of 
the nearest-neighbor Heisenberg kagome-lattice antifer- 
romagnet, we argue that the spectacularly strong quan¬ 
tum effect of the quasiparticle breakdown in an almost 
classical system should be applicable to a variety of other 
fiat-band frustrated spin systemsWe also remark on 
a useful q-dependence of the dynamical structure factor, 
which is characteristic to the non-Bravais lattices and al¬ 
lows to select spectral contributions of specific branches 
in the portions of the q-space. 

The paper is organized as follows. In Sec.[n|we present 
a consideration of the harmonic theory for several com¬ 
mon anisotropic models on the kagome lattice, explicate 
details of the diagonalization procedure, and show results 
of a representative calculation within the harmonic ap¬ 
proximation. In Sec. HI, anharmonic terms of the models 
are derived and the results of the ground-state selection 
calculations are presented. The spectral properties of the 
kagome-lattice antiferromagnets are also given a detailed 
exposition. Technical aspects of the derivation of the 
quartic terms are given in Appendix \A\ 


II. LINEAR SPIN-WAVE THEORY 

To set the stage, we provide the linear spin-wave con¬ 
sideration of the isotropic nearest-neighbor antiferromag¬ 
netic Heisenberg model on the kagome lattice following 
the approach of Ref. [5l We continue with various ex¬ 
tensions of the model, which are either relevant to real 
materials or allow to explore the role of quantum effects 
in a wider parameter space. These extensions include 



(b) f ^ 

t t 

^ t 

t t 

V3xV3 I ^ 


FIG. 1: (Color online) (a) q = 0 and (b) \/3 xy^ spin config¬ 
urations. In (a), primitive vectors of the kagome lattice and 
numbering of sites within the unit cell are shown. 


the anisotropic XXZ model, models with the single-ion 
and the out-of-plane Dzyaloshinskii-Moria anisotropies, 
and additional further-neighbor exchange terms. Subse¬ 
quently, the important steps of the diagonalization proce¬ 
dure that will be essential for the non-linear terms consid¬ 
ered in Sec. [TIT] are exposed. In the end of this Section, 
results of the calculations of the on-site magnetization 
within the linear spin-wave theory for the XXZ model 
are presented as an example. 


A. Nearest-neighbor Heisenberg model 

Within the spin-wave treatment, nearest-neighbor 
Heisenberg antiferromagnet on the kagome lattice with 
the Hamiltonian 

n = j^Si-Sj, ( 1 ) 

{ij) 

is assumed to be in an ordered state with spins forming 
a coplanar 120° structure. Summation is over bonds and 
are the sites of the lattice. Aligning the z-axis on 
each site along the direction of the ordered moment and 
directing the ^-axes out of the ordering plane transforms 
the Hamiltonian 0 to a local spin basis 

^ = JY .+ CO®) 

{ij) 

+ sin% {StS^ - SfS^)) = Si ® S,- , (2) 

(b') 

where Oij = Oi — Oj is an angle between two neighbor¬ 
ing spins and we have introduced “matrix” product G) of 
spins with the matrix 

( cos Oij 0 — sin Oij \ 

0 10 (3) 

sin^^j 0 QOS Oij I 

as a shorthand notation. 

We choose the unit cell of the kagome lattice as an 
up-triangle, see Fig. and the primitive vectors of the 
corresponding triangular Bravais lattice as 

<5i = (l,0), <52= 53 = <52-<5 i, (4) 
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all in units of 2a where a is the interatomic distance. 
The atomic coordinates within the unit cell are pi = 0, 
P 2 = 5^3, and pa = -^Si. 

Then, changing summation over bonds to summation 
over unit cells and the atomic index, i with 

a = 1-3 enumerating atoms within the unit cell, the 
Hamiltonian (§ becomes 

n = J^Si.,®(S 2 ,, + S2/-3) (5) 

i 

(8) + Ss^^+i) + S2/ (8) (S3/ + S 3 ^^+ 2 ) 5 

where the product (8) Sa',£' is according to ([^, and 
i ± n = ± with the coordinate of the unit cell 

= miSi + m2^2- 


with the matrix 

/ 0 C 3 Cl \ 

Ak = C 3 0 C 2 I , ( 10 ) 

\ci C 2 0 y 

and shorthand notations c^ = cos(g'^) with qn = 'k - Sn/2. 
One can rewrite this Hamiltonian as 

^2 = 5I^kHkXk-3J5, (11) 

k>0 

with the vector operator 

= (^l,k’ ^Ik’ ^3,k’ ^ 1 ,-k’ ^ 2 ,-k’ %,-k) 


B. Harmonic spin-wave approximation 


It should be noted that although a coplanar state with 
spins on each triangle in a 120° structure minimizes the 
classical energy of 0, such a state is not unique and 
the manifold of them is extensive .1^ However, it is also 
clear from the Hamiltonian in the local basis 0 that the 
linear spin-wave theory is the same for any state from 
this degenerate manifold,® because cos Oij = —1/2 for 
any pair of spins in such state and terms do not 

contribute to the harmonic order of the 1/S expansion. 

Thus, we introduce Holstein-Primakoff representation 
for spin operators 




a,£ 


S V 2-5 - (6) 


into 0 and, keeping only quadratic terms, obtain a har¬ 
monic Hamiltonian for three species of bosons 


7^2 — ^ (^ 2 ,£^ 2 ,£ ^ 3 ,£^ 3 ,£ 

£ ^ 

^1,£{^2,£ ^ 2 ,^- 3 ) + ^1,£{^3,£ %/+l) 

+ + ®3/+2) +h.C. 


1 

+ 8L' 


Vn 


“a,I 


^*krQ, 


(7) 


+ ^2/-3) ^1,£{^3,£ + %/+l) 

+ “2,<'(®3,£+ “ 3 ,^+ 2 ) +h.C. I . 


Performing the Fourier transformation according to 

1 


(8) 


and the 6x6 matrix Hk 


where 

Ak = i+iAk, Bk = ^Ak, (14) 


and i is the identity matrix. 

For a moment, we will confine ourselves to the eigen¬ 
value problem of 7 / 2 - Because of an obvious commuta¬ 
tivity of the matrices Ak and Bk, the eigenvalues of Hk 
are straightforwardly related to their eigenvalues, and, 
in turn, are determined by the eigenvalues A^y^k of the 
matrix Ak, so that the spin-wave excitation energies are 

£,,k = = ^JScOuM , (15) 


with the frequencies cj^y^k = y/(l — Ajy^k/2)(1 + A^y^k) and 
1 3 

Azy,k = 1 ^ ^ • (1^) 

Thus, the problem of the diagonalization of 7^2 in ^ is 
reduced to the eigenvalue problem of Ak in (10). From 
the characteristic equation for Ak one finds 


|Ak - A| = (A + 1) (A2 - A - 27k) = 0 


(17) 


where 7 k = C 1 C 2 C 3 is introduced and factorization is per¬ 
formed with the help of a useful identity 

Cl + C2 + C3 = 1 -h 2C1C2C3 , 


where rQ ,7 = R^ + pQ, and N is the number of unit cells, 
we obtain the linear spin-wave theory Hamiltonian 


772 = 2J5' 



^a/3 


4 ^ 


^a,k^/3,k 


3 
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K^{‘ 


'a,k 


a 


t 

/3 -k 



(9) 


which holds once the cosine arguments satisfy q 2 = Qi + 
qs. Thus, the A^y^k eigenvalues are 

= “1 5 "^2(3),k ~ 2 ^ \/l + 87 k^ . (18) 

Of the resultant spin-wave excitations one is completely 
dispersionless and has zero energy, referred to as the “flat 
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mode,” and two are “normal,” i.e. dispersive modes, 
which are degenerate in the Heisenberg limit 

^l,k = 0, ^2(3),k = — 7k • (19) 

The nat ure o f the flat mode has been discussed 
previously!!!^ Generally, such modes owe their ori¬ 
gin to both the topological structure of the underlying 
lattices that facilitate spin frustration and the insuffi¬ 
cient constraint on the manifold of spin configurations. 
Physically, they correspond to the localized, alternat¬ 
ing out-o f-pl ane fluctuations of spins around elementary 
hexagonswhich do not experience a restoring force in 
the harmonic order in the Heisenberg limit, hence their 
energy is zero. In the following, various anisotropies lift 
the energy of such a flat mode, but preserve its flatness. 


XXZ model 


Next, we consider an extension of the Heisenberg 
model on the kagome lattice to the XXZ model with 
anisotropy of the easy-plane type, 0 < A < 1. The origi¬ 
nal motivation for this extension, see Ref. [9l was that the 
degeneracy among the 120° coplanar states in this model 
remains the same as in the Heisenberg model. This has 
allowed us to extend the parameter space and to study 
the effect of quantum fluctuations in the ground-state se¬ 
lection without lifting d egen eracy of the classical ground- 


state manifold, see Sec. HI for more detail.l^ 


In the case of the XXZ model, the plane for the copla¬ 
nar 120° structure is chosen by the anisotropy. In the lo¬ 
cal spin basis of Q with ^-axis out of the ordering plane, 
the XXZ addition to the Heisenberg model reads as 


(ij) 


( 20 ) 


A now straightforward spin-wave algebra of (20) leaves 
the structure of the harmonic Hamiltonian in intact, 
yielding corrections to the Hamiltonian matrix in (13) 


( 21 ) 


<5Ak = <5Bk = AA Ak . 


The spin-wave energies for the XXZ model in ( |15| ) are 
now with = y^(1 — Ajy^k/2)(1 + AA^^^k) and give 


ei,k = 2JS^/3(l^A)/2, (22) 

for the flat mode, which is now at a finite energy, and 

£2(3),k = 2JS^ 1 - A7k - (1 - A)(1 ± x/l + 87k )/4. 
for the dispersive modes. 


D. Single-ion anisotropy 


Instead of the XXZ correction (20), an alternative way 
of generating easy-plane anisotropy is to add a positive 


single-ion term 


6H = Dj2iS!y 


(23) 


where y is the out-of-plane axis in the local spin basis of 
(© as before. This term, in a complete similarity to the 
XXZ case, gives zero contribution to the classical en¬ 
ergy, does not affect cubic anharmonicity, and does not 
contribute to the degeneracy lifting of the 120° mani¬ 
fold through the quartic terms .1^ Its contribution to the 
harmonic Hamiltonian (13) is also simple 


(5Ak = 5Bu = -i 


(24) 


where d = D/J and I is the identity matrix. Thus, once 
again, the eigenvalue problem of 7^2 reduces to the eigen¬ 
value problem of Ak, resulting in the spin -wave frequen¬ 
cies cj^y^k = y/(l — Ajy^k/2) (1 + d + Ajy^k), ^ud the spin- 
wave energy of the flat mode 


ei,k = 2JS^M/2, 

while energies of the dispersive modes are 


(25) 


£2(3).k = 2JS^J 1 + d - 7k - c;(l ± y 1 + 87k )/4, (26) 

which should be compared to the XXZ results above. 

There is a high degree of similarity of the single-ion 
anisotropy model ( [M| ) and its results to the XXZ case, 
the most important being no degeneracy lifting among 
the 120° manifold of classical states at the harmonic level 
of approximation. However, there is also an important 
difference. If analyzed in real space in the local basis 
the structure of the spin-flip terms is different in the two 
models. In particular, the single-ion term (23) creates 


spin flips that are purely local while the spin-flip hop¬ 
ping and other amplitudes responsible for degeneracy¬ 
lifting remain independent of anisotropy D. Therefore, 
from the point of view of the real-space perturbation the¬ 
ory, described in Ref. [H the minimal order in which an 
effective degeneracy-liftiM interaction is generated is dif¬ 
ferent in the two models.^ We will address this difference 
in Section in more detail. 


E. Dzyaloshinskii-Moria interaction 

Important correction to the Heisenberg spin model 0 
that commonly occurs in real magnetic materials with 
the kagome structure is the anisotropic Dzyaloshinskii- 
Moriya (DM) interaction. 

The out-of-plane DM term, see Fig.|^ is the main per¬ 
turbation to the Heisenberg Hamiltonian of the S' = 5/2 
kagome-lattice antiferromagnet Fe-jarositej ^^Fd g =1/2 
materials herbertsmithit^!^ and vesignieite,!!^ and other 
systems .1^ This anisotropy differs from the ones con¬ 
sidered above in a number of aspects, but allows for a 
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FIG. 2: (Color online) Directions of the out-of-plane DM vec¬ 
tors. Arrows on the bonds show the ordering of the Si and 
Sj operators in the vector-product in (27). 


very similar analytical treatment, mainly because of its 
nearest-neighbor nature. 

The antisymmetric Dzyaloshinskii-Moriya interaction 
is generally written as 


^'Hdm = ^ Dij • (Si X Sj). 

{'ij) 


(27) 


In order to determine the DM vectors Dij one has to 
specify the order of sites in the vector product, which 
may be represented by the bond direction from the first 
to the second spin in each pair. A convenient choice 
consists of ordering spins uniformly along the chains, see 
Fig-i Then, the symmetry analysis yields a unique pat¬ 
tern of the DM vectors o rthogonal to the kagome plane 
as shown in Fig. Note, that the in-plane compo¬ 

nents of the DM vectors are strictly forbidden once the 
the kagome plane coincides with a mirror crystal plane.^^ 
Apparent alternation of the DM vectors between up 
triangles with = ( 0 , 0 , —D^) and down triangles with 
'Dij = ( 0 , 0 , Dz) is partly fictitious, because it is a conse¬ 
quence of the chosen bond gauge, in which the two types 
of triangles are circumvented oppositely, Fig.[^ The DM 
interactions on the two triangles favor the same sense of 
spin rotation or chirality. Hence, for a given sign of Dz^ 
the DM term (27) selects one of the two q = 0 struc¬ 


tures with positive {Dz > 0 ) or negative {Dz < 0 ) chi¬ 
ralities, yielding energy gain = —V3\Dz\S‘^ per site. 
On the other hand, for the x state contributions 
from up and down triangles come with opposite sign and 
cancel out. In the following, we assume Dz>0 and, con¬ 
sequently, perform the spin-wave expansion around the 
q = 0 state with positive chiralities, see Fig. It is also 
straightforward to see that the DM term does not in¬ 
duce additional canting and preserves the 120 ° magnetic 
structure in the x-y plane. 

Consider the DM term (27) on the (1, 2 ) bond in Fig.|^ 
written in the rotating local spin basis ([^ 

^ D,smei2iSISl+SIS^) 

+D, cos 012 - S^St) , (28) 
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FIG. 3: (Golor online) Energies of the three magnon modes 
for the J-Dz spin model with the out-of-plane DM interaction 
Dz/ J = 0.06 along a representative path in the Brillouin zone. 


where Oij =0i — Oj , and the first term contributes to the 
classical, harmonic, and quartic orders of the 1/S expan¬ 
sion, while the second contributes in the cubic order. 

Since the DM term concerns only the nearest-neighbor 
pairs of spins and because in the q = 0 state all DM 
bonds contribute identically, the overall structure of the 
harmonic part of the Hamiltonian remains the same as in 
Then, some algebra yields the harmonic Hamiltonian 


7^2 in the form (13) with the DM contributions 


(5Ak = ( I - ^ Ak 


<5Bk=^Ak, (29) 


where (Im = V^Dz/J and Ak is unchanged from ( 10 ), 
again reducing the eigenvalue problem of the harmonic 
spin-wave theory to the one of Ak, already solved in 
(17) and ( p^ . Then, the spin-wave frequencies for the 
problem with the out-of-plane DM interaction (p^ are 

= a/(1 + du) (1 ~ A^y,k/2) (1 + djM + With 

Ajy^k from (p^, the fiat mode energy is 


£i,k = 2JS^Mm (1 + dM) /2 , (30) 

and the energies of the dispersive modes are 

^2(3),k = 2 JS'y/r^T'dM (31) 

X ^1 + (iM - 7k - (l ± \A + 87 k ) /4 , 
which are similar to the XXZ 


and the single-ion 
(25), (26) results and are in agreement with Ref.| 2 T| The 
energies of the magnon modes 5jy^k are illustrated in Fig.[^ 
for the DM t erm D; 2 /J = 0.06, a value relevant to the Fe- 
jarosite.l^^^ 

Opposite to the previous extensions, the DM term con¬ 
tributes to the cubic anharmonicities, see (28). Using the 


second term in (28) for one bond and extending it to the 
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entire lattice we find for the q = 0 state 


% = ^ ^ sin% (aja]aj + h.c.) 


(32) 




whose structure is identical to the cubic term from the J- 


part of the Hamiltonian (47) considered in Sec.|III[ Thus, 


the out-of-plane DM interaction in the q = 0 state simply 
renormalizes cubic vertices by a factor (1 + 


F. Small-J 2 expansion 


In the Heisenberg kagome-lattice antiferromagnet, ad¬ 
ditional next-nearest-neighbor coupling J 2 lifts the de¬ 
generacy of the 120 ° manifold of classical states and se¬ 
lects between the q = 0 and \/3 x \/3 ground states 
It also introduces a dispersion into the “flat mode” en¬ 
ergy and thus was invoked to reproduce experimentally 
observed k-dependence of such mode in Fe-jarosite!^ We 
note, however, that quantum fluctuations beyond the 
harmonic order also generate effective J 2 interactions , 1 ^ 
and thus could be the source of the same k-dependence. 
As we show in Sec. in the dispersion of the flat mode is 
particularly important for the quasi-resonant spin-wave 
decays. Below, we consider the effect of small J 2 pertur- 
batively. Other types of small interactions can be taken 
into account in a similar fashion. 

We first point out that the network of the second- 
neighbor bonds forms three independent kagome lattices. 
Second, these bonds connect only spins from different 
sublattices of the original lattice, see Fig. There¬ 
fore, contribution of the J 2 term to the harmonic spin- 
wave Hamiltonian has the same structure as the nearest- 
neighbor Hamiltonian 


5% 
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t 

a,k 


a 


t 

/3,-k 



(33) 


where instead of Ak the matrix is 

/ 0 4 4 \ 

^k = I ^3 0 ^2 I 5 (34) 

V 4 4 0 J 

and we use the shorthand notations c[ = cos(g 3 + g 2 ), 
C 2 =cos(g 3 -gi), C 3 = cos(gi+g 2 ) with qn = k-Sn/2. 

Therefore, the diagonalization of the harmonic part of 
the J — J 2 Hamiltonian requires diagonalization of the 
matrix Ak = Ak + where j 2 = J 2 lJ- Generally, this 

can be done numerically, the approach taken in Ref. [ 21 ] 
in the analysis of the Fe-jarosite spectrum, with the an¬ 
alytical results given only for the high-symmetry points. 
However, having in mind the problem of spin-wave exci¬ 
tation in large-S' kagome-lattice antiferromagnet, for the 



FIG. 4: (Color online) Evolution of the “flat mode” upon 
switching on additional interactions with parameters shown 
in the legend. Solid lines: energies of the “flat mode” for the 
J — J 2 Heisenberg (lower curve) and out-of-plane DM models 
(upper curve), using approximate expressions in Eqs. (37) and 
( |39| ). Dashed lines are the energies for the same models at 
J 2 = 0. Dotted lines (virtually indistinguishable from solid 
lines) are exact results, without using the expansion in J 2 . 


physical range of interest J 2 ^ J one can make analytical 
progress using an expansion in jf 2 . 

Because the most important qualitative effect of J 2 
is the dispersion of the flat band, we will ignore small 
corrections to the “normal” modes. Expanding charac¬ 
teristic equation for the matrix Ak in j 2 we obtain 

(A + 1) (A^ — a — 27k) = 2jf2 (A/i(k) + / 2 (k)), (35) 
where /i(k) = c^ci + 402 + C 3 C 3 , 

/ 2 (k) = 4 C 2 C 3 -h C 2 C 1 C 3 -h C 3 C 1 C 2 . 


Then, the “fiat mode” Ai = —1 eigenvalue is modified as 

^,. = A.+i.AS‘>, ASyAfflyffi). (36) 


Corrections to A 2 and A 3 can be obtained similarly. 


Having this perturbative correction to Ai in (36) allows 


us to obtain the fiat mode energies in various models. We 
list some of the answers below. 

(i) J — J 2 Heisenberg model: 


£i.k = 2J5^3 j2(1 + a7)/2 + 0{jy^). (37) 


(ii) J — J 2 XXZ model. Here we assume that the 
anisotropy A is the same in both exchanges: 


ei.k = 275-^3/2+j2(l-A^^72) (38) 

X \Jl - A + j 2 (l + A^^^) + 0 (j|). 
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(iii) J — J 2 out-of-plane DM model. 


£i,k = 2JSyJz{l + dM)/2 + J 2 (l - a2/2) (39) 

X \ldn + j2(l + -h 0(j|), 

where (Im = V^Dz/J as before. This is the model which 
was used to describe the spectrum of Fe-jarosite.^^ Our 
results agree exactly with the expressions for the high- 
symmetry points provided in Ref. [211 The advantage of 
our approach is that it is fully analytical in the entire 
Brillouin zone. 

Effects of the second-neighbor exchange J 2 and the 
DM coupling Dz on the dispersion of the “flat mode” 
^i,k are summarized in Fig. for representative values 
J 2 /J = 0.03 and DzjJ = 0.06 tha t are motivated by the 
experimental data for Fe-jarosite p^ * ^^ * the J 2 = 0 ener¬ 
gies are also shown. In the same figure the energies ob¬ 
tained via numerical diagonalization of the matrix Ak are 
shown. The corresponding results are indistinguishable 
from the approximate results of Eqs. (37) and (39) on 
the scale of our plot. 


G. Two-step diagonalization 

For the above three cases, the XX the single-ion, 
and the out-of-plane DM models, the harmonic spin-wave 
theory includes diagonalization of the same matrix Ak- 
Here we elaborate on details of this general procedure 
and provide the formalism, which is essential for treating 
the non-linear terms in all these models and is identical 
for all considered cases. 

Following Ref. [5j diagonalization of Ak implies a two- 
step diagonalization procedure of I-L 2 ia the form (13). 
Its eigenvectors wj = (re, 


[;^,i(k),n;^, 2 (k),n;^, 3 (k)) obey 
AkW^ = A^,kW^ (40) 

and can be found explicitl}®^ 


w^,(k) = — 

Tv 


C1C2 + AjyCS 
C1C3 + AjyC2 


(41) 


The final step is the conventional Bogolyubov transfor¬ 
mation for each individual species of the d-bosons 




Kk- 

-Kk 

— 

2 

_ 1 / 

'4 


2\ 

. UJi 


^zy,k ^ 


- 1 , = 


B, 


i/,k 


^zy,k 


(43) 


(44) 


where the eigenvalues cj^y^k, ^zy,k: and were obtained 
for each of the models in previous sections. The impor¬ 
tance of this two-step procedure will be apparent in the 
discussions of the non-linear terms in Sec. imi 


H. Ordered magnetic moment 

Here we discuss the dependence of the ordered mag¬ 
netic moment on anisotropy parameter and on the value 
of the spin S. While we only consider the XXZ model,!^ 
the results are expected to be similar for the other 
anisotropies considered above. It should be noted that 
this calculation can only qualitatively estimate the sta¬ 
bility of the Neel order, because it only includes the 
“diagonal” quantum fluctuation for a given state, com¬ 
pletely neglecting the “off-diagonal” tunneling processes 
between different states within the manifold. On the 
other hand, such processes should be exponentially sup¬ 
pressed for larger spins 

Within the linear spin-wave theory, magnetic moment 
on a site that belongs to the sublattice a is reduced by 
zero-point fluctuations: {S)a = S — (a^ ^). Convert- 

and Bo- 


ing aa to du and then to using unitary 


golyubov (43) transformations one arrives to 




Vk • 


(45) 


/x,k 


Since all three sublattices are equivalent, symmetrization 
of (45) gives 




37V ^ ’ 

^t,k 


(46) 


with = \J{c\C 2 -h A^C 3 ) 2 -|- (A^ - ciy+ (C 1 C 3 -I- A^C 2 )^. 

These eigenvectors define a unitary transformation of 
the original Holstein-Primakoff bosons to the new ones 


^a,k — 


= E 


w, 


(k) = E (k) a„.k > (42) 


such that the harmonic Hamiltonian 'Ho in the form of 
^ with Ak’s and Bk’s from (14) with ( [2l| ), (24), or (29) 
is turned into three independent Hamiltonians 


H2 = 2JS 




k^zy,k 






«k<-k + h.c.) 


zy,k 


with from (44). Calculations of the magnetization 
M = (S) /S and tEe {S) = 0 Neel order boundary in the 
S—A plane are performed taking the 2 D integrals in (46) 
numerically using ( [4^ for the XXZ model. 

Quantum suppression of the ordered moment vs 
anisotropy A is shown in Fig.|^a) for two values of spin. 
Linear spin-wave theory suggests a disordered state near 
the Heisenberg limit for all spin values because quantum 
correction diverges for A —1 due to vanishing energy 
of the “flat mode.” The critical value 1 —Ac ~ 0.047 for 
5' = 1/2 is compared in Fig. [^a) with the result for the 
DM coupling Dz^c = ^-^ (rescali ng 1—A c VSDz.c)^ found 
by exact diagonalization (ED).^i^^ Since the DM term 























FIG. 5: (Color online) (a) Magnetization M—{S)/S vs A for 
S—1/2 and S — ^/2. Square is the ED result for 5'=l/2 with 
the DM interaction, (b) Intensity plot of the magnitude of 
M. Solid line is the {S) =0 Neel order phase boundary in the 
S'—A plane on the log-log scale, dashed line is the asymptotic 
approximation for it, see text. 


spin-wave expansion to conclude the formal development 
of the theory. Then the cubic vertices for q = 0 and 
Xy/S states allow us to proceed with calculating order- 
by-disorder fluctuating corrections to their ground-state 
energies for the XXZ and single-ion anisotropy models. 
Both models demonstrate a quantum phase transition 
between the q = 0 and the x y/3 states as a func¬ 
tion of anisotropy parameter. Hence, both cases present 
rare examples of the quantum order-by-disorder favoring 
a different state from the one selected by thermal fluctua¬ 
tions, the latter c hoosing the v^Xy/3 structure regardless 
of the anisotropy.l^l^l^^®^lI^K^ 

We then proceed with a calculation of the decay- 
induced effects in the structure-factor 5(q, cj) within the 
DM model with J 2 . While this calculation is aimed at 
giving a detailed account of the spectral properties of a 
specific kagome-lattice antiferromagnet described by this 
model, S = b/2 Fe-jarosite, the outlined scenario should 
be applicable to a wide variety of the other flat-band frus¬ 
trated spin systems.l^^M^ We also note that the structure 
factor q-dependence allows to “filter out” spectral con¬ 
tributions of specific modes in the portions of the q-space 
while highlighting the others; a useful phenomenon char¬ 
acteristic to the non-Bravais lattices. 


A. Cubic terms 


suppresses tunneling processes within the manifold, it is 
reasonable to compare ED with spin-wave theory results 
to evaluate the accuracy of the Neel order boundary. One 
can see a qualitative agreement with ED and a quanti¬ 
tative exaggeration of the extent of the ordered phase, 
expected for the spin-wave approach. 

Near the Heisenberg limit, one can neglect non- 
divergent terms in ( [4^ and find an asymptotic expression 
for the Neel order boundary from {S) = 0 ^ S' —Ai/bcJi, 
where Ai = 3/4, see (16), and cji = ^3(1 — A)/ 2 , see 
( [^ , leading to 1 — Ac ~ 1 /96S^, which is shown in 
Fig. [^b) together with the result of the numerical in¬ 
tegration in (46) and demonstrates an exceedingly close 
agreement with it. 


III. NON-LINEAR SPIN-WAVE THEORY 

In Sec.|n|we have considered three anisotropic models 
of the kagome-lattice antiferromagnets in the harmonic 
approximation and outlined the approach to taking into 
account other terms, such as further-neighbor exchanges, 
perturbatively. In this Section, we derive the nonlinear, 
cubic and quartic terms of the spin-wave 1 /S expansion 
and then exemplify their role in the ground-state selec¬ 
tion and in the spectral properties of the kagome-lattice 
antiferromagnets using representative cases. 

Below, we first obtain cubic and quartic terms of the 


In the XXZ, single-ion anisotropy, and Heisenberg 
models, (§ with or without the anisotropies in ( [2Q| ) and 
(23), the terms that lead to the cubic anharmonic cou¬ 
pling of the spin waves are identical and originate from 
the SfS| part of ([^. These terms are also the only 
ones that are able to distinguish between different 120 ° 
spin configurations in these models by virtue of contain¬ 
ing sin Oij = ±^/3/2 for the clockwise or counterclockwise 
spin rotation. In the bosonic representation they yield 

'H3 = j \j^ Z. (4 


where Oij = Oi — Oj is the angle between two neighboring 
spins as before. Clearly, the spin-wave interaction result¬ 
ing from this term has different amplitude for different 
spin patterns. For the DM model, the DM term (27) also 


contributes to the cubic anharmonicity, but its structure 
for the q = 0 state is identical to (47), see (32), so it 


gives a simple change of the overall factor in the vertex 
J ^ J -^Dz / \/3- 

Using ( [T7[ ), we now detail the derivation of the cubic 
vertices for the main contenders for the ground state, the 
q = 0 and the \/3 x \/3 states. We begin with the q = 0 
pattern, for which 7/3 in (UtI can be rewritten as 


% = - jy ^ Z cos(g; 3 „)a^_qa^ k«/ 3 .P + h- c., 


Q:/3,k,q 


( 48 ) 

























9 


where is the Levi-Civita antisymmetric tensor, p = 

k + q, Q/3a — Q ■ P/3a-) p^a ~ P /3 Po 


The unitary transformation (42) in (48) yields 


'35' 


% = - ^ E E KQ + h.c., (49) 

k,q yfir] 

where p = k + q and the amplitude is 

-^qkp = E COs(5/3a) W^,a{ci)u!f,,i3{k)Wr,,0{p). (50) 

a/3 


Finally, the Bogolyubov transformation (43) gives 
1 1 




3! //V 
1 1 


EE^qTp 4q4,k4,-p+ll-C- (51) 


k,q u/J-r] 


E E Kk-p ^tq4,k^^.p+h-c., (52) 

k,q u/J-r] 

with the vertices for the “source” and the “decay” terms 

IS = -•'V? • ®Xp = -•'V? *5’p ■ (53) 

where the symmetrized dimensionless vertices are 


~ ^qkp('^i^q ^ '^z^q)('^Aik'^?7p + '^/xk'^?7p) 

+ l^Cq(“A*k + ^'Mk)( '^zyp'^r 7 q T '^zyp'^? 7 q) 
4“ -^pqk(^i?P 4” '^77p) ('^z^q'^/xk 4“ '^zyq'^/xk) ; 

and 


(54) 


^qk^p ~ ^qkp('^i^q + '^zxq) ('^t^k'^? 7 p + '^^tk'^? 7 p) 

4“ -^i^pqC'^/xk 4“ '^^k) ('^zyp'^?7q 4“ '^zyp'^?7q) (55) 

4“ -^pqkC'^TVP 4“ '^77p) ('^zyq'^^k 4“ '^zyq'^/ik) 5 

where we have used the symmetry property . 

Repeating the same calculation for the ^/3 x ^/3 state 
we obtain identical expressions for the cubic spin-wave 
Hamiltonian and corresponding vertices, but with differ¬ 
ent amplitudes F^^^ expressed as 

-^qkp = * E sin(g;3„) W^^a{ci)Wf,,0{k)Wr,,0{p). (56) 

a/3 

Explicit expressions for the unitary transformation eigen¬ 
vectors of the matrix Ak and of the parame¬ 

ters of the Bogolyubov transformation are instrumental 
in deriving analytic expressions of the cubic anharmonic 
terms. We also note that for all models considered here, 
the functional form of the cubic vertices (54) and (55) is 


identical, with all the differences hidden in the expres¬ 
sions of the Bogolyubov parameters and 'r’^k from 
(|44|). 


The role of the cubic terms in the ground-state selec¬ 
tion and in the spectrum of the kagome-lattice antiferro- 
magnets is discussed below in Sec. |III C] and Sec. HID 


B. Quartic terms 

In the Holstein-Primakoff bosonic representation of 
spin models, quartic terms originate from the 5f5J, 
5f 5J, and 5f 5| parts of the Hamiltonian. In the models 
considered here, quartic terms do not help to differen¬ 
tiate between different states of the 120° manifold, but 
lead to the overall ground-state energy shift and to the 
Hartree-Fock corrections to the spin-wave energies. Be¬ 
cause of the coplanar 120° spin configuration, the formal 
expressions for these contribution show a close similarity 
to the ones for the triangular-lattice Heisenberg model, 
see Ref. n Therefore, we simply list the quartic parts of 
the Hamiltonians of the XXZ model and the DM term 
together with the expressions for the ground-state energy 
shift in the former model and for the spin-wave energy 
correction for the latter model. Technical details are of¬ 
fered in Appendix [A| 

The quartic terms in the XXZ model, ([|) and 


are 


^4 = 4 ^ + nj)aiaj + h.c.) (57) 

{ij) ^ 


2A - 1 


(aj(ni + nj)ai + h.c.) — riin^ , 


where rii = ata^ and we omitted the a indices of the 
bosonic variables for brevity. 


The quartic contribution from the DM term (27) is 

^4, 


DM 


^ ^ ((^i + nj)aiaj + h.c.) (58) 

{ij) ^ 

-h ^ {a]{ni + nj)ai + h.c.) - . 

By means of the Hartree-Fock decouplin^^ outlined in 
Appendix one can obtain contribution of the quartic 
terms to the 1/5 series expansion of the ground-state 
energy of any non-collinear spin model 

F = Fci + , (59) 

where the first term is the classical energy of order 0{S^), 
the second is the harmonic correction from 
and the last two are from the nonlinear cubic and quartic 
terms, both 0(1). While we discuss the cubic part of this 
expression later, the classical and harmonic energy terms 
(per spin) of the series in (59) for the XXZ model are 


Ec\ 




I 5 


^t,k 


(60) 


with the frequencies from ([22|). The quartic member of 

(61) 


the series (|^ for the same model is 

A 2 


5Ea = -J 


■ m 


-(2A + 1)( 


nA ■ 


mS\ 

~Y) 


(2A — 1) (^nm -\- 


A(5\1 

) 
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where the Hartree-Fock averages are introduced in a stan¬ 
dard manner 

n = (ata^), m = (ata^), A = 5 = (a^), (62) 

and are given in Appendix \A\ 

In the same order of expansion, quartic terms lead 
to the Hartree-Fock corrections to the linear spin-wave 
Hamiltonian 

sn, =2JJ2 ^ (4,k4,-k + h. c.), 

iy,k 

which yields the Hartree-Fock part of the 1/S contribu¬ 
tion to the spin-wave energy 



^(4) _ O 7 


^v,k 


(63) 


Here we give expressions for and for the 

Heisenberg model (§ with the DM anisotropy ( [2^ as 
an example 

5^.,k = Cl + = - 2 C 4 - (64) 


FIG. 6: (Color online) Difference of the ground-state energies 
of the q = 0 and \/3 X\/3 states per spin. Upper inset: energy 
correction for the q = 0 (squares) and \/3 x \/3 state 

(diamonds). Dashed vertical lines mark the transition. Upper 
data points/lines in the figure and inset are for the single-ion 
anisotropy model, § with ( |23| ), and the lower are for the 
XXZ model, ([^ with ( [2o| . Lower inset: diagram for 
term (65) in the energy expansion (59). 


where Ci = Ci — dMDi and the constants are the linear 
combinations of the binary Hartree-Fock averages in ( [6^ 
and can be found in Appendix We note that since 
the “flat-mode” eigenvalue Ai^k = —1, the quartic 1/S 
correction to its energy, in (63), is also necessarily 
momentum-independent in all rnoaels considered here. 
Therefore, it is a contribution of the cubic terms which 
is going to introduce a fluctuation-induced dispersion in 
the flat mode in the same order of the 1/S expansion. 


C. Ground state selection 


Here, we discuss the role of cubic terms (51) in the 


ground-state selection in the XXZ and the single-ion 
anisotropy models. As was mentioned above, in the case 
of these models, <§ with ( [2Q| ) and ( [M| ), classical and 
harmonic terms are not able to lift the degeneracy in the 
manifold of coplanar 120° structures. It is also easy to 
see from <§ that the quartic terms are also unable to 
differentiate between 120° states, leaving the cubic term 
as a sole source of the quantum order-by-disorder effect 
to this order in 1/S. The same is also true for the single¬ 
ion term in ( [2^ . 

The second-order energy correction from the cubic 
terms (51) resulting in 6E^^^ in ( [H^ is represented by 
the diagram in the lower inset of Fig. and is given by 


5E^^^ = 


EE- 


KmU- 




(65) 


where the energy is per spin and N is the number of 
unit cells. Summation over magnon branches gives 


twenty seven individual contributions of which ten are 
distinct. With the formal expression for the source ver¬ 
tices in (54), the energy correction in (65) is identical for 
the XXZ and the single-ion cases, with the differences 
in the expressions for the spin-wave energies in (22) 

and (25) and ([26|, and in the u^k and v^k parameters 
(44) in vertices (|54p for the corresponding models. 

Using an explicit form for the cubic vertices for the 
q = 0 and \/3 x \/3 spin states, we performed a high- 
accuracy numerical integration in Eq. ( [^ and studied 
the quantum order-by-disorder effect. Results of these 
calculations are presented in Fig. The quantum se¬ 
lection of the q = 0 state over the \/3 x \/3 counterpart 
for the large planar anisotropy (1 — Ac) [Dc/J] >0.3 was 
highlighted in our previous work on the XXZ model.l^ 
This qualitative conclusion is in contrast with the usual 
behavior, in which quantum fluctuations lead to the same 
ground state that is selected by the thermal fluctuations. 
Indeed, for the classical kagome-lattice antiferromagnets 
in both Heisenberg and XY limits, thermal fluctuations 
select the \/3x\/3 magnetic structure as the leading 
instabiht}^^^^^^^^®^ contrary to the behavior of the quan¬ 
tum model. 

Here we complement the above result by the analy¬ 
sis of the single-ion anisotropy model (23). While the 
overall trend in particular, 

the transition point between the states in both models 
in Fig. I^are very similar, there is a quantitative differ¬ 
ence. As we have noted previously, the local real-space 
structure of the degeneracy-lifting terms is different in 
the XXZ and single-ion models, implying a higher-order 
real-space path needed for generatiM the ground-state 
selection in the latter model.^ Fig. [^provides a strong 
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A A 


FIG. 7: (Color online) Energy per site of the q = 0 state in 
the XXZ model. Dotted line is for the classical term E n\ = 
— JS^, dashed is Ed + SE^‘^^ in and solid line is for (59) 


with 6E^^^ and 6E^^^ from ( [dlj ) and ( [6^ . Square is the DMRG 
energy for the spin-liquid state.® (aj^ =1/2, (b) S = 1. 


support to this point, as the energy difference is smaller 
for the single-ion model by a factor 2 — 8 for the range of 
0.3 < (1 — A)[D/J] <1, in a qualitative agreement with 
the parameter of the real-space expansion being ^ Ijz 
(z = 4, coordination number).!^ 

We conclude the discussion of the ground state with 
Fig.[^ which shows the energy of the q = 0 state in ( [^ vs 
A for the XXZ model for two representative spin values 
S'= 1/2 and S=l. Classical and harmonic contributions 
to the ground-state energy are also indicated by dotted 
and dashed lines. Note that the energy difference shown 
in Fig.j^would be nearly invisible on the scale of the plot 
in Fig, p We also note that although 8E^^^ and SE^"^^ 
from ml) and ( [6^ represent the entire contribution of 
the 0(1) order in the 1/S expansion of the ground-state 
energy (59), their divergences at A ^ 1 do not cancel 
completely, thus signifying a breakdown of the standard 
1/S expansion at the Heisenberg limit because of the 
vanishing energy of the flat mode.^ 


D. Spectrum and decays 


We now turn to the spectral properties of the kagome- 
lattice antiferromagnets. The goal of our consideration 
is twofold. First, we would like to highlight an unusual 
spectral property that has to be present in a wide va¬ 
riety of frustrated spin systems with excitations featur¬ 
ing flat branches.® Second, we give a detailed account 
of such spectral p roper ties in a specific model that de¬ 
scribes Fe-jarosite.®® Because of that we concentrate 
on the out-of-plane DM-anisotropy model, (|^ with ( [^ , 
of the S'= 5/2 nearest-neighbor kagome-lattice antiferro- 
magnet, which closely describes Fe-jarosite in the range 
of small Dz. One can expect the results of this considera¬ 
tion to be similar to the ones for the XXZ and single-ion 


anisotropy models given the similarity of their harmonic 
spectra and anharmonic terms, even though the ground- 
state selection is more subtle in the later models. 


Formalism and a qualitative discussion 


Regardless of the model, using standard diagrammatic 
rules with the cubic terms in (51) and (52) produces the 


spin-wave self-energy in the form 






q,i /?7 ^»7.k-q 

IC-k-q.kP 


-iS 


( 66 ) 


^ + ^zy,q + ^r7,k—q — iS j 


where the first and the second terms are the decay and 
the source self-energies. Because of the summation over 
the magnon branches u and r] in the decay and source 
loops, there are nine terms in the sum in (66), only six 


of which are distinct. Note that for the DM model one 
has to change J^ J-\-Dz/V^ in the vertices. 

Taken on-shell, uj = the self-energy (66) repre¬ 


sents a strictly 1/S correction to the magnon energy from 
the cubic terms, 5]^,k(^/x,k) =0{S^). The other term of 
the same order in the 1/^'-expansion is from the quar- 
tic terms, which was obtained for the DM model in 


(63). Then, the magnon Green’s function for the branch 
/i can be written as 


^(k, w) = w - - eG “ , 


(67) 


which, generally, allows to evaluate the spectral function 
A^(k,ci;) = — (l/7r)ImG^(k, cj) of the corresponding spin- 
wave branch ji. 

Since we are interested in the large-S' limit and in 
the resonance-like decay phenomenon in the spectrum, 
the following simplification can be used. Given the off- 
resonance character of the source term, the frequency- 
independence of the quartic terms, and the large-S limit 
of the problem, one can neglect the real part of the 1/S 
corrections to the spectrum as a first step, and approx¬ 
imate the self-energy in (66) by its on-shell imaginary 
part, i.e. S^(k,w)RizImS^,e^,k) = -*r^,k, with 


^ E I^Xq;kP<^ (^M.k - e..q - £,.k-q), (68) 


q,iy?7 


where the summation is over magnon branches of the 
decay products. Using the dimensionless vertices in (55) 
and frequencies in (30) and (31), one can rewrite (68) as 


^t,k 

IT 


37^ I q,k-q;t 


^ (^^t,k ^zy,q ^?7,k—q) (69) 


which is explicitly independent of the spin value S. As is 
mentioned above, the summation over magnon branches 
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in (69) contains six distinct terms, or potential “decay 


channels.” A particular decay channel may or may not be 
contributing to the decays, depending on the kinematic 
conditions discussed below. 

Before we divert our attention to a specific model, the 
following qualitative consideration can be made. The 
form in ( [ 6 ^ suggests that T^^k = which is small 

compared to = 0{JS) in the large-5' limit. Thus, 
generally, one can expect a regular, perturbative 1/5- 
effect of broadening of the higher-energy part of the 
dispersive branches due to decays into the lower-energy 
states.^ However, a spectacular exception to this rule 
must occur if both of the decay products, branches u 
and 7 ^, belong to flat modes with uji = const. Clearly, 
this decay channel produces an essential singularity in 
T^^k and S^(k,cj) at the energy equal to twice the flat 
mode energy, the effect that can be seen as a resonance¬ 
like decay. Therefore, a special care must be exercised 
in this case as any residual dispersion of the flat mode 
becomes crucial in regularizing such a singularity. 

In fact, the way of removing this singularity is offered 
by the same 1/5 fluctuations due to cubic terms, as the 
self-energy ( 66 ) also yields a dispersion for the flat mode. 


cation of the problem is in having three different magnon 
branches, which modifies these condition to 


This can also be interpreted as a fluctuation-generated 
further-neighbor J 2 spin-spin inter actions,!^ and such 
a dispersion of the fiat mode has been observed in Fe- 
jarosite.l^ Still, this scenario implies that the fluctuation- 
induced flat-mode bandwidth, Wi oc 0 ( 1 ), is of the order 
1/5 compared to the bandwidths of the normal disper¬ 
sive modes, 1 F 2 ( 3 ) oc 0(5). 

Introducing this 1/5 dispersion for the fiat modes 
in (69), suggests that the near-resonance decay-induced 


broadening of the dispersive mode within the energy 
window of the width 21 Fi near 2 co’i should now scale 
as T^^k = 0(J5). This is the same energy scale as 
the spin-wave dispersion itself without any obvious ad¬ 
ditional smallness. Therefore, our analysis implies a 
very strong broadening, likely eliminating spectral weight 
nearly completely from the respective energy range even 
when spin 5 is large, providing an example of a spectac¬ 
ular quantum effect in an almost classical system. 

Thus, in theory, as 5 ^ oc, we predict that an anoma¬ 
lous broadening and a wipe-out of the spectral weight 
should remain in the spectra of the flat-band frustrated 
systems, albeit in the energy window of order 0{J) while 
the spectrum width grows as 0(J5). In practice, we 
argue that such a spectacularly strong quantum effect, 
leading to the quasiparticle breakdown with characteris¬ 
tic termination points and ranges of energies dominated 
by broad continua, must be present in an almost classical 
5 = 5/2 kagome-lattice Fe-jarosite. 


2. Decay channels in Fe-jarosite 

An extensive analysis of the kinematic decay condi¬ 
tions in representative models of t he fru strated spin sys¬ 
tems has been provided previously) ^^ * ^^ * Here, the modifi¬ 


^/i,k — T ^?7,k—q i 


(70) 


SO that every decay channel {z^, 77 } has to be analyzed 
separately. That is, for each branch /i up to six channels 
can be contributing independently to the decay rate 


r^,k — ^ rk,p 


>{u,r]} • 


(71) 


While the decay conditions are model-specific, Fe- 
jarosite offers a commonplace scenario: a predomi¬ 
nantly nearest-neighbor Heisenberg system with sublead¬ 
ing symmetry-breaking anisotropic term, which is re¬ 
sponsible for lifting the flat mode to a finite energy. We, 
therefore, analyze the decay channels for this represen¬ 
tative situation for a small Dz = 0.06J, used in fitting 
spin-wave spectrum of Fe-jarosite.l^ The harmonic dis¬ 
persions for this value of Dz are shown in Fig. [^and in 
the inset of Fig. [^a), and we use them in our analysis. 

We use the same numeration of the branches as be¬ 
fore: l=fiat mode, 2 =gapless mode (Goldstone branch), 
and 3=gapped dispersing mode. While the flat mode 
does have one active decay channel into two Goldstone 
modes, 1 ^ { 2 , 2 }, for a small range of momenta near 
the center of the Brillouin zone, k < kc where kc is the 
intersect point of the branches, see Fig. [^a), its effect 
is truly minor and the main interest is in the decays of 
the dispersive modes. From the picture of harmonic en¬ 
ergies, it is easy to see that the energy conservation in 
( [^ can be satisfied for three (four) decay channels out 
of the six for the mode 2(3), one being the “resonance¬ 
like,” 2(3) ^{1,1}, mentioned above, and two (three) are 
“regular.” The latter are 2(3) ^ {1, 2} and 2(3) ^ {1, 3}, 
with an additional channel for the mode 3^{ 2 , 2 }. 


Using (69), we calculate contributions of the individual 
decay channels, Fk^^^jzy^Ty}, shown in Figs.l^a) and (b), 
for modes 2 and 3 and for a representative k- cut of the 
Brillouin zone. One can see that for 2(3) ^ {1,1} h 
gives a (5-peak at 2co’i, that 2(3) ^ { 1 , 2 } channel has a 
threshold at uJi and 2(3) ^ {1,3} at 2uJi. In Fig. [^a), 
the resonance decay rate is scaled down by 10 “^ and 
would otherwise dwarf the rest even for the small artificial 


broadening of the ^-function in (69). 


A somewhat unexpected finding is a strong suppres¬ 
sion of the 3 ^ { 2 , 2 } and the 3 ^ {1,1} channels. A 
closer analysis have identified different origins for them. 
For the 3 ^{2, 2} channel, decay products involve long- 
wavelength excitations from the Goldstone branch, which 
provides an extra smallness in the decay amplitude com¬ 
pared with the “regular” ones. For the 3^{1,1} channel 
the situation is more subtle. One can demonstrate that 
the vertex ^q^k-q k carries an extra Dz compared with 
the vertex for the 2 ^{ 1 , 1 } channel, yielding a factor of 
about 1/100 in the decay rate. The physical reason for 
this suppression can be hypothesized. Quasiclassically, 
the mode 2 is the “in-plane” mode, while modes 3 and 
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FIG. 8: (Color online) Individual contributions to the decay 
rate in (69) from the decay channels in units of J, 

(a) for mode 2, (b) for mode 3 in for the DM model. Dz — 
0.06J, k is along the FY direction, and artihcial broadening 
(5 = 0.002 has been used, k- values that correspond to ui and 
2(jUi are indicated by the arrows. Results for the “resonant” 
decay channel 2^{1,1} are scaled down by the factor 10~^. 
Inset: harmonic frequencies Ct;a,k along the representative cuts 
of the BZ for all three modes from Fig. 


1 are “out-of-plane.” Then the in-plane mode can cou¬ 
ple naturally to the two out-of-plane modes, while the 
out-of-plane mode has hard time splitting into two. 

Thus, the suppression of the resonant decays of the 
mode 3 seems to be due to a subtle cancellation in the 
structure of the corresponding vertex. We note that for 
a realistic case of Fe-jarosite, other symmetry-breaking 
terms are also present, so one can expect the subtle can¬ 
cellations of more symmetric models to be violated. 


3. Spectrum of Fe-jarosite 

We now finally approach the realistic description of 
the Fe-jarosite spectrum. We note that more details of 
this discussion is offered elsewhere.lHl As was mentioned 
above, the essential singularity in the dispersive modes 
is naturally removed by the residual dispersion of the 
flat mode. The experimentally observed dispersion of 
the flat mod^^ has been fitPH by introducing small next- 
nearest neighbor exchange J 2 = 0.03J, ignoring its pos- 



FIG. 9: (Golor online) Lower curve with the shading is the on- 
shell F 2 ,k from (69). Dashed line is the linear spin-wave theory 
energy of the gapless dispersive mode, £: 2 ,k, from (31). Shaded 
area shows the half-width boundaries of a lorentzian peak, 
e: 2 ,k±F 2 ,k- Dotted and solid lines are different approximations 
for the renormalization of the real part of the self-energy in 
(66), see text. Parameters are as shown in the plot. 


sible quantum origin. Since we do not perform a fully 
self-consistent calculation here, the same approach suf¬ 
fices for the removal of the singularity in the decay rate 
( [6^ . We, thus, modify the fiat-mode dispersion accord¬ 
ing to (39) and ignore other corrections from the J 2 term. 

In Sec. HIE we elaborate on the possible way of sep¬ 
arating contributions of different modes in the neutron¬ 
scattering structure factor, which allows us to concen¬ 
trate on individual modes. Therefore, our Fig. sum¬ 
marizes the effects of decays on the gapless dispersive 
mode only, as they are most pronounced in it. We have 
used parameters shown in the figure with the value of 
J = 3.34meV from the previous work.^ 

The lower curve shows the on-shell F 2 ,k obtained from 
(69) with the fiat-mode dispersion from (39), and it in¬ 
cludes all three decay channels discussed above. The 
dashed line is the linear spin-wave energy, 52 ,k, with the 
shaded area around it representing 52 ,k i r 2 ,k 5 i-e., half¬ 
width at the half-maximum boundaries of a lorentzian 
peak. We have also taken into account renormalization 
of the real part of the self-energy in (66), with the dotted 
line showing the 1/S on-shell result and the solid line in¬ 
cluding effect of self-consistency by taking into account 
imaginary part from F 2 ,k in calculation of ReE. One 
can see that the resultant effects of the nonlinear terms 
on the real part of the spectrum are relatively minor, in 
agreement with the discussion after Eq. (67). 

We complement Eig. by an intensity map of the 
spectral function, A 2 (k, cj), for the same dispersive mode 
along the same representative path in the Brillouin zone 
and for the same parameters, see Eig. Dashed and 
dotted lines are the spin-wave results from Eig. I^and are 
guides to the eye. The magnon self-energy is approxi- 
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FIG. 10: (Color online) Intensity plot of the spectral function 
A 2 (k,(jj) of the gapless mode in units of {2SJ)~^ along the 
same path and for the same parameters as in Fig.[^ lines are 
from the same hgure, see text. A small broadening S/2SJ = 
0.002 has been added to F 2 ,k in (68). Left panel inset: decay 
diagram. 


mated by its on-shell imaginary part, as discussed above. 
The upper cut-off of the intensity of the spectral func¬ 
tion corresponds to the maximal height of the peaks in 
the non-resonant decay region in Fig. and translates 
into the broadening Fk ~ 0.73meV for the Fe-jarosite, 
easily resolvable by the modern neutron-scattering ex¬ 
periments. 

Given the large spin value, S' = 5/2, we estimate that 
the ordered moment in Fe-jarosite should be nearly 90% 
of its classical value, see Fig. It is then very natural to 
assume that the spectral properties should be fully de- 
scribable by the harmonic spin-wave picture, the point of 
view taken in Ref. [21] Indeed, our Figs. and [^demon¬ 
strate that the spin-wave excitation is sharp below the 
flat-band energy £i^k and acquires only an infinitesimal 
width for the energies above it. However, there is a sharp 
threshold at twice the bottom of the flat band minimum, 
265^Jf, above which there is a very strong broadening, 
reaching about one third of the bandwidth, si gnif ying an 
overdamped spectrum.El Above 2s^^^ in Fig. 10 there is 
a broad energy band with the features that fcok like a 
rip in the spectrum, consistent with the missing spectral 
weight in the experimental data. This threshold singu¬ 
larity is also remarkably similar to the spectral signatures 
of the quasiparticle breakdown pheno meno n in quantum 
Bose liquids and S = 1/2 spin-liquids 

Although at the energies above twice the top of the 
flat-band maximum, 2e^^, the decays are “regular,” 
i.e., occurring due to other, non-resonant channels, they 
are still providing a measurable broadening to the spec¬ 
trum, so it is instructive to compare the two regions. 
The “regular” decays result in the maximal values of F 
of order 0.2 — 0.3 J, an agreement with the similar ef¬ 
fects in triangular-lattice antiferromagnet^I^ and other 



FIG. 11: (Golor online) Intensity map of A 2 (k,ct;) in units of 
(2SJ)~^ vs k throughout the Brillouin zone for for the same 
set of parameters as in Fig. I^for uj = 18meV. The upper cut¬ 
off of the intensity scale corresponds to the broadening F 2 ,k of 
0.73 meV for the Fe-jarosite values of S and J. Dashed black 
lines are peak positions from the linear spin-wave theory. 


frustrated spin systems.!^ The maximal broadening in 
the resonant-decay region is F « 1.7J for the Fe-jarosite 
model, which is larger than the effect of the “regular” 
decays by a factor close to 5 (= 2S). This is in a remark¬ 
ably close accord with the qualitative argument on the 
scaling of the resonance decay rate with the spin value, 
provided after Eq. (69) above. 


While a more detailed description of the spectral fea¬ 
tures of Fe-jarosite is given elsewherewe would like to 
highlight here a different, perhaps more dramatic view on 
the drastic transformations in its spectrum in the range 
of energies that falls within the resonant-decay band. In 
Fig. we show an intensity plot of the constant-energy 
cut of the spectral function, A 2 (k,co’), of the same dis¬ 
persive mode as in Figs. and [To] for the energy 18meV, 
which is close to the middle of the resonant-decay region. 
The dashed lines show the expectations from the linear 
spin-wave theory of where the contours of sharp, well- 
defined peaks should have occurred. Instead, one can 
observe strong deviation from such expectations, char¬ 
acterized by a broadening, massive redistribution of the 
spectral weight into different regions of k-space, together 
with the multitude of intriguing features, which are re¬ 
lated to the van Hove singularities in the two-particle 
density of states of the flat-band decay products. [^Alto¬ 
gether, Figs. [9 , jT , and IT offer a convincing evidence of 
the highly non-trivial and very strong quantum effects in 
the the dynamical response of a nearly classical flat-band 
kagome-lattice antiferromagnet, which are facilitated by 
the nonlinear couplings. 
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FIG. 12: (Color online) Kinematic formfactors: for the 

(a) flat, (b) dispersive gapless, (c) dispersive gapped mode; 
Fiq for the (d) flat, (e) dispersive gapless, (f) dispersive 
gapped mode. in (e) diverges at some F points as 


E. Dynamical structure factor 


Inelastic neutron scattering cross-section is directly re¬ 
lated to the diagonal components of the dynamical struc¬ 
ture factor, or the spin-spin dynamical correlation func¬ 
tion, given by 

df 

5*°*°(q, co) = J - {t)S\), (72) 

where io refers to the laboratory frame {xo,^o,^o} and 


qio — qio 


(73) 


involves summation over the spins a in the unit cell. 

Because of the coplanar spin configuration in the con¬ 
sidered kagome-lattice antiferromagnets, transformation 
from the laboratory reference frame of io lo the local 
spin basis of ([^ yields a mix of different diagonal and 
off-diagonal terms in the structure factor,!^ which conve¬ 
niently separate into the in-plane and out-of-plane parts 
of uj). Assuming equal contribution of all three io 

components to the neutron-scattering cross sectioiP^ and 
using the mapping of spins on bosons ^ allows one to 
perform a straightforward 1/S ranking of different con¬ 
tributions to the structure factor, in which the transverse 
components are, as usual, dominate the longitudinal and 
mixed terms. 

The subsequent algebra involves the two-step transfor¬ 
mation (42) and (43) from the Holstein-Primakoff bosons 


to the quasiparticles, yielding the leading contributions 
to the structure factor as directly related to the spin-wave 
spectral functions Ajy(q, cj) 


5in(out)(q^^) « ^ (74) 

V 



FIG. 13: (Golor online) A 2D map of the inverse lifetime 
from (77) in units of J in the q-space. A small artificial 
gaussian q-broadening with a = 0.027r/a was used to mimic 
instrumental resolution. The scale is given in the inset and 
the maximal value of F§^q^ 1.7 J corresponds to about 6meV 

for Fe-jarosite. 


where we introduced kinematic formfactors 

+ Vuqf (1 - Ruq) , 

= f Kq - (1 + 2i?,q) , (75) 

with 

Rtyci = ^ 51 ^<^,a(q)^!^,a'(q) • (76) 

a^a' 


An important property of the kinematic formfactors 
in (75) is their q-dependence. They are modulated in 
the q-space and are typically suppressed in one of the 
Brillouin zones while are maximal in the others. This 
effect is characteristic to the neutron-scattering in the 
non-Bravais lattices and is, in a way, similar to the effect 
known as the Bragg peak extinction for the elastic scat¬ 
tering in such lattices. Because of that property, one may 
be able to focus on a specific excitation branch without 
intermixing contributions from the others by selecting a 
particular component of the structure factor in a partic¬ 
ular Brihon in z one. 

Our Fig. shows 


for three branches of ex¬ 
citations in the Fe-jarosite model. It demonstrates that 
such a “filtering out” can be quite useful. For exam¬ 
ple, the out-of plane component of the structure factor, 
cj), should be dominated in one of the three dis¬ 
tinct Brillouin zones by the spectral function A 2 (q, cj) of 
only one of the dispersive modes. This feature can be 
utilized in the neutron-scattering experiments. 

Lastly, we highlight another quantitative way of ana¬ 
lyzing 5(q, cj). Assuming that one can focus on a par¬ 
ticular excitation branch as mentioned above, one can 
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suggest that a faithful representation of the inverse life¬ 
time (linewidth) F^q from (68) and its distribution in 
q-space across the Brillouin zone can be obtained from 
the moments of the dynamical structure factor aP^ 


Tq = - (w5(q,w))2, (77) 

where (cc;’^5(q, cj)) = J uj'^S{q^,uj)duj are the moments of 
the structure factor and normalization (5(q, u;)) = 1 is 
assumed. Then, this experimentally unbiased procedure 
would allow extracting a 2D q-map of the quasiparti¬ 
cle broadening. We demonstrate the effectiveness of this 
approach in Fig. which shows an example of such a 
map derived using the procedure in from a gaus- 
sian form 742(q, cj), i.e. the gaussian function with a 
maximum at 52,q and the width r 2 ,q. As expected, the 
extracted map of Fq is nearly identical to the map of 
F 2 ,q itself. However, for a more natural lorentzian form 
of A 2 (q, cj) one can show that the extracted map corre¬ 
sponds to Fq oc >/ci;niaxr2,q, whcrc cjmax Is the Upper limit 
of the integration over uj in the moments in ( [77| ). Impor¬ 
tantly, this result is still providing a direct information 
on the quasiparticle broadening map, albeit on a differ¬ 
ent scale. Therefore, aside from demonstrating at which 
momenta the decays are most intense, the suggested pro¬ 
cedure also provides another way of “fingerprinting” of 
the broadening due to resonance-like decay into the flat 
modes. 


IV. CONCLUSIONS 

In summary, we have provided a systematic considera¬ 
tion of the nonlinear 1/S expansion of several anisotropic 
models of the kagome-lattice antiferromagnets. We have 
demonstrated the role of the nonlinear terms in the quan¬ 
tum order-by-disorder selection of the ground state and 
presented a strong evidence of the rare case of quan¬ 
tum and thermal fluctuations favoring different ground 
states in two of these models. We have provided a de¬ 
tailed analysis of the excitation spectrum of the S' = 5/2 
iron-jarosite to illustrate our proposed general scenario 
of drastic transformations in the spectra of the flat-band 
frustrated magnets. Our study calls for further neutron¬ 
scattering experiments in these systems. 
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Appendix A: Hartree-Fock corrections 


The quartic terms in (57) yield a correction to the 
ground-state energy given by the four-boson averages, 
which are decoupled into the products of the binary 
Hartree-Fock averages (62) using Wick’s theorem 


{a\aia}-aj) = , (AI) 

(ata^a^aj) = 2nA + md, (atata^ai) = 2nm + Ad, 


where the Hartree-Fock averages are obtained similarly 

to the calculation of the staggered magnetic moment in 

(45). Namely, for the on-site averages n = 

and 6 = following the transformations (42) 

and (43) from to and to and using equivalence 

of all three sublattices one arrives to the symmetrized 

expressions 


= (4,^ Ydk> 

I 


^ ^ Y (A2) 




with and from (44). For the nearest-neighbor 
two-site averages, the same transformations lead to 


^ = («b«/3,r) = V Y //3«(k)dk > 


^t,k 




(A3) 


^t,k 


where // 3 a(k) = cos{kga) a P, 

are the components of the eigenvector in the 


transformation ( [4^ , and = kpg^ with = Pg—Pa 
as before. Although it seems that the two-site averages 
may depend on the choice of a and /3, one can verify that 
all three possible combinations of a ^ pairs yield the 
same answer. 

Quartic terms also yield Hartree-Fock contributions 
to the linear spin-wave dispersions via (63). The cor¬ 


responding constants for the Heisenberg and DM terms 
in (64) are 


3A m 3d 

Ci = -n+ — - C 2 = -m+ — 

^ ^ 3n 5 ^ 3m A 

C3 = -A+--j. = 



£>3 = A- 


n 

2 


m 

Y’ 

d 

4’ 


D2=m-- 


n 

2 ’ 



2 ’ 

(A4) 
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